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Abstract. Very recently [Phys. Rev. E 82, 021921 (2010)] a simple mechanism 
was presented by which a molecule subjected to forced oscillations, out of thermal 
equilibrium, can maintain quantum entanglement between two of its quantum degrees 
of freedom. Crucially, entanglement can be maintained even in the presence of very 
intense noise, so intense that no entanglement is possible when the forced oscillations 
cease. This mechanism may allow for the presence of non-trivial quantum entanglement 
in biological systems. Here we significantly enlarge the study of this model. In 
particular, we show that the persistent generation of dynamic entanglement is not 
restricted to the bosonic heat bath model, but it can also be observed in other 
decoherence models, e.g. the spin gas model, and in non-Markovian scenarios. We 
also show how conformational changes can be used by an elementary machine to 
generate entanglement even in unfavorable conditions. In biological systems, similar 
mechanisms could be exploited by more complex molecular machines or motors. 
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1. Introduction 

The last couple of years have witnessed the birth and explosive development of a new 
research area, namely quantum biology. For many years the possibility of non-trivial 
quantum effects in biological systems was thought to be impossible due to the large 
amount of noise that characterize such systems. During the last few years however, very 
surprisingly, a host of experimental and theoretical works have pointed to the possible 
presence of such effects. Notable results concern the existence of quantum coherence 
in photosynthesis PQElElllll5lEllZllSll3llIDlinilI2llI3!aiid quantum effects in the so 
called "chemical compass" [H [151 ESI EH US CES] . 

Of particular interest is the possible existence of non-trivial entanglement in 
biological systems. The reason is that entanglement has dramatic effects, in particular 
it leads to a vastly increased capacity of information processing. Indeed, entanglement 
is at the core of quantum computation. The capacity of quantum systems to perform 
many other tasks is also significantly increased by entanglement. Given these benefits 
it is very tempting to suspect that biological systems may make use of entanglement. 
On the other hand, the existence of entanglement in biological systems was thought to 
be even more unlikely than that of mere quantum coherence because entanglement is 
extremely fragile and almost impossible to survive in the presence of noise. Very recently, 
however, it was realized [20] that the intuition about the fragility of entanglement 
was based on considering the behavior of quantum systems at thermal equilibrium. 
Biological systems, however, are not at equilibrium - they are open driven systems, 
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far from thermal equilibrium. Such systems are capable of "error correction" and may 
actively maintain entanglement in presence of noise. This raises the possibility that 
non-trivial entanglement may actually be present in biological systems. 

In [2TJ a simple mechanism was presented by which a molecule subjected to forced 
oscillations, out of thermal equilibrium, can maintain quantum entanglement between 
two of its quantum degrees of freedom. Crucially, entanglement can be maintained even 
in the presence of very intense noise, so intense that no entanglement is possible when the 
forced oscillations cease. A potential biological system in which such an effect may take 
place is that of complex molecules (such as proteins) undergoing forced conformational 
changes. At present, no such system has been experimentally observed, for the simple 
reason that no experiment has yet been carried out to look for it. To be able to do that 
it is imperative to first gain a better theoretical understanding of the model. This is 
the scope of the present paper. 

The paper is organized as follows: In section [2] we present the abstract model 
of a moving "two-spin molecule" with time-dependent coherent Hamiltonian; then we 
investigate the system behavior and the entanglement dynamics under the influence of 
two distinct environments: In section [3] we consider a bosonic heat bath representative 
of a large number of harmonic oscillators, and study its influence both in the Markovian 
and non-Markovian regime. To relax the working assumption of periodic oscillations, 
we further introduce stochasticity in the classical motion of the molecule, and include 
an additional source of noise to illustrate the robustness of the entanglement generation. 
In section [4] we extend our study to the spin gas model. The final section [5] is devoted 
to some concluding remarks. 

2. How conformational molecular motion can generate entanglement 

In the following, we consider bio-molecular systems that undergo controlled 
conformational changes, using a semi-quantal picture. We assume that the 
conformational motion of the molecular structure can be described, to a very good 
approximation, classically. In addition, however, we identify selected quantum degrees 
of freedom, embedded in this "backbone structure" and localized at several classically 
defined sites, which interact with a strength dependent on their distance (see Fig. pi). 




Figure 1. (color online) Conformational changes of a bio-molecule [22], induced e.g. 
by the interaction with some other chemical, can lead to a time-dependent interaction 
between different sites (blue circles) of the molecule. See also j20]. 



In the following, we consider two quantum degrees of freedom which are forced to 
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move along some predetermined trajectories and interact via dipolar coupling. Each 
of these degrees of freedom shall be represented by a quantum mechanical two-level 
system. For easy reference, we call such a model system, a "two-spin molecule" [21]. It 
can be regarded as an elementary example of a molecular machine which, as we shall 
see, is able to persistently generate entanglement. As such, it also presents a building 
block which could be exploited in more complex bio-molecular motors. 



.!'(/) X 2 (l) 




Figure 2. A two-spin bio- molecule undergoes conformational changes such as 
harmonic oscillation or stochastic motion. Both the spin-spin interaction strength 
J and the background field B are position dependent. (Fig. taken from [2T] ) 

In more detail, the mechanism is described as follows: The conformational changes 
in the molecule force the two localized spins (or, more generally, two-level systems) 
to come close or move apart, inducing an effective time-dependent coupling, both 
between the two spins and between each spin and a local field, as described in Fig. 
Furthermore, it will be assumed (this is essential) that the de-coherent environment 
is so noisy that no static entanglement between the two spins can exist in the static 
steady state, which is achieved when the molecular configuration is fixed to any of the 
possible ones. As we shall see in the following, if the molecular configuration oscillates 
in a periodic way, cyclic generation of fresh entanglement can persist even in the long 
time limit. This feature will persist even in the presence of an additional source of 
dephasing, which would generally eliminate all coherence from the system. Finally, we 
shall generalize these considerations to more realistic scenarios involving non-Markovian 
environments and stochasticity in the classical motion. 

In the following, whenever we refer to spins, we may imagine spin-1/2 magnetic 
moments interacting via a dipolar coupling and subjected to a background magnetic 
field. To perform our study, we further assume that the two spins are coupled with the 
usual Ising type interaction which depends only on the distance d between the spins, and 
that there exists a local magnetic field which induces a (Zeeman) energy level splitting. 
The energy splitting is also assumed to depend only on the inter-spin distance d, which 
thus plays the role of an effective parameter for the overall molecular configuration. 
Explicitly, the Hamiltonian for the system is 

H M (d) = J(d)a^aW + B {d){aU + <r< 2 >), (1) 

where ai a \ a* are Pauli operators of the a-th spin, oo (d) = 2B(d) is the local level 
splitting, and we express everything in units of the Planck constant (U = 1). Since 
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the molecule undergoes conformational changes, the distance d between the two spins 
becomes a function of time: d = d(t). The explicit time dependence of Hj^id) = 
HM(d(t)) is determined by the specific choice of the molecular motion: Initially, we 
consider forced periodic oscillations, while we later include an additional stochasticity. 
With this model, the molecular motion leads to an effective time dependent Hamiltonian 
for the system that we can write as 

H M (d(t)) = H M (t) = J(t)a^ + B(t)(a^ + a?\ (2) 

where, here and in the following, we use the short notation J(t) = J(d(ij) for all the 
quantities which depend on the distance d and so implicitly on the time t. For notational 
simplicity, we introduce the energy scale 



S(t) = a/ AB 2 (t) + J*(t) (3) 

in terms of which the eigenenergies of the system Hamiltonian H^it) and the 
corresponding eigenstates can be written as 
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(4) 



with 

n m ^ = ^v_pv . (5) 

It can be seen that the energy spectra of the system Hamiltonian HM(t), and in 
particular the ground state |e ), depends on both the interaction strength J it) and 
the local field B(t), and thus is also time dependent. In case of weak spin-spin 
interaction J(t) B(t), e.g. when two spins are spatially separated from each 
other, the system Hamiltonian can be approximated as Hm ~ B(t)(ai 1 ' 1 + ai 2 ^), and 
the ground state approaches the separable state |11). As the spins move closer and 
the spin-spin interaction increases, the ground state changes according to |e (t)) in 
Eq.Q. In particular, the entanglement contained in the ground state |e (i)) is a 
monotonically increasing function of the ratio j^L and can be quantified by the 
measure of concurrence [23]. Concurrence is zero for separable states and reach its 
maximum value C = 1 for maximally entangled states; here it can be calculated as 
C(|e (t)}) = | (eo(t)\<r v ® a y |e *(t)) | = 2^/(1 + V 2 ) = [1 + 4(H) 2 ]" 1 , where |e *(t)> is 
the complex conjugate of state |eo(t)) (expressed in the computational product basis, 
already used in Eq.Q ). In the meantime, the relative energy gap between the ground 
state and the first excited state (ej — eo)/ (£3 — eo) decreases with the ratio J{t)/ B{t). 

Before taking the environment into account, it is worth looking at the adiabatic 
coherent evolution of the oscillating molecule, which will facilitate our discussions in the 
following sections. If the system Hamiltonian changes slowly in time, the instantaneous 
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eigenspaces of HM{t) will undergo an independent evolution. The widely used condition 
that assures the adiabatic dynamics for closed quantum systems is 

,MM|,<< 1 , kjti. 

For a detailed discussion concerning the validity of the adiabatic approximation we refer 
to [2HI2S]- For the time-dependent system Hamiltonian introduced in Eq.(|2]) one easily 

gets 



|ei(*)> = |ca(*)> = (7) 



and 



\eo(t)) = h(t)) , \e 3 (t)) = T ^_ \e (t)). (8) 



< 2E(t). (9) 



f^|6 3 (t)), |^(f)) = -2- 

1 + r] z 1 + if 

Thus the adiabatic condition can be explicitly expressed as 

I V 

I I + r] 2 

Whenever the adiabatic condition is satisfied, the molecule that is initially in an 
eigenstate of the system Hamiltonian will remain in the corresponding eigenstate, that 
is to say 

|e*(0)) U(t,0) |e fe (0)) ~ e -</<N*M* e **<*) \e k {t)} , (10) 
where f Q ek(s)ds is the dynamical phase and 9 k {t) = i (efe(s)| e'k(s))ds represents a 



geometric phase. One can check from Eq.([^|8]) that (ek(t) \ e&(i)) = is fulfilled Vfc and 
all times. Thus no geometric phase is accumulated, i.e. 9k(t) = for all k, and the 
adiabatic evolution can be expressed as follows 

|e fc (0))^e-^o^(^| efe (t)). (11) 



3. Bosonic thermal bath 



The coherent part of the spin dynamics due to the time-dependent system Hamiltonian 
has to be incorporated into a more complex description which takes into account 
the presence of an environment. The actual environment of biological systems is 
very complex and probably sophisticated, and far from being thoroughly understood. 
However, one can adopt certain representative and well-investigated decoherence models 
to gain first insights into the open dynamics of bio-molecules. We first demonstrate 
the essential idea by considering the familiar decoherence model of spins coupled with 
independent thermal baths consisting of a large number of harmonic oscillators. The 
Hamiltonian of the individual thermal bath coupled to the a-th spin (a = 1, 2) is 

k 
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where a k a and the creation and annihilation operators for the fc-th bath 

oscillator mode and the corresponding mode frequency. The coupling between the 
spin and the thermal bath is described by the following interaction Hamiltonian 

HSB = £ H$> = E ^ E ^.«( fl M + °U (13) 

where is the coupling amplitude to mode k. 

Within the Born-Markov approximation, and by assuming that the adiabatic 
condition for closed quantum systems Eq.Q is satisfied, we have derived a time 
dependent Lindblad-type quantum master equation [2T] for the system dynamics under 
the rotating wave approximation, which reads as 

j t p{t) = -i[H M (t), p(t)} + V(t)p(t) = C(t)p(t), (14) 
with the following dissipator 

V(t)p(t) = E r -W {A a Ht)}p(t)AUto(t)] - Ai[u(t)}A a [u(t)}p(t)} + h.c. 

Ct,U)(t) 

A a [uj{t)\ = fcWkf x !*(*)> • ( 15 ) 

Aij(t)=oj(t) 

In other words, the time dependence of the system Hamiltonian introduces an implicit 
time dependence into the Lindblad generators A a [u;(£)], as well as into the rates r w ( t ) 
(a sort of one-side Fourier transform of the thermal bath correlation function, see [2l]). 
Notice that, as already discussed for Eq.Q, the time dependence of H M (and so also 
T> and C) is actually due to its dependence on the distance d between the spins, from 
which the short notation, e.g., £(t) = £(d(t)). 

The bosonic thermal bath is characterized by its spectral density function J {to) 
which contains information on both the density of modes in the frequency domain and 
the strength of their coupling to the system. In the following, we choose an Ohmic 
spectral density J {to) = ^ue^ w ^ c , with cutoff frequency uj c . The memory time of the 
thermal bath strongly depends on the cutoff frequency and so does its Markovian or 
non-Markovian behavior [26]. In particular, if u c ^> u, where u is the energy scale of 
the system transition frequencies, the reservoir loses its memory quickly, and is called 
Markovian. The corresponding correlation function exhibits a sharply peaked kernel in 
the limit u c — > oo. In this case, one can simplify the spectral density as J{u) = ku with 
k = 4-- We thus can write the Fourier transformation of the thermal bath correlation 
function as = Ku(t)(l+N u a),p) where N^uyp = l/(e w w^ — 1) represents the bosonic 
distribution function at inverse temperature /3. 

In order to establish a comparison between the static and dynamic cases, we 
first consider the situation in which the two-spin bio-molecule is fixed in a static 
configuration. Since the two spins are coupled to independent thermal baths at the same 
temperature, the corresponding quantum master equation mixes the energy eigenstates 
in the sense that the molecule is driven towards its thermal equilibrium state [2S] at the 
reservoir temperature 1//3. For each spin-spin interaction strength J(d) = J and local 
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magnetic fields B(d) = B, the entanglement of the corresponding thermal equilibrium 
state, defined by C(d)p th = 0, is quantified by [21] 

C(jH h ) = -= max {0, — ^ sinh - cosh (J/3)}, (16) 
Z 1 + r] z 

with 2 = 2cosh(£/3) + 2 cosh (J/3) the partition function. After some calculations, it 
can be seen that the first derivative of C with respect to /3 is always non negative 

dC(p th )/d/3 > 0, (17) 

which means that the static thermal entanglement decreases as the temperature 
increases; in particular, static entanglement can not survive if the temperature is higher 
than a certain critical temperature 1//3 C - In other words, if the molecule stays in a fixed 
configuration, entanglement will finally disappear for high enough temperatures, e.g. 
see Fig. (|3]). In this sense, the bosonic heat bath, even though it may not be a very 
realistic model of a biological environment, plays a useful role in our investigation. It is 
an example of an unfavorable environment, which cannot introduce any entanglement 
by itself. 




Figure 3. Static entanglement in the thermal equilibrium state at different inverse 
temperatures /? for B(x) = B(0) — x and J(x) = J(Q) + x with B(0) — 1.2 and 
J(0) = 0.1. Here f3,B,J can be varied independently, but it is their relative value 
which determines the dynamical behavior. Therefore, here and in the following, we 
use dimensionless quantities obtained by setting the thermal energy k^T as the unit 
of the energy scale (so that j3 — 1) with respect to which all the other parameters are 
determined, for example B = 1 indicates that HB/k^T = 1. 

Here we address the key question: What happens if one takes into account the 
molecular motion? and in particular: Can entanglement possibly build-up through the 
classical motions of the molecule? The answer is affirmative and we demonstrate in 
the following that entanglement can indeed persistently recur in an oscillating molecule, 
even if the environment is so hot that any static thermal state is separable, i.e. the 
temperature T > max{T c } for all possible molecular configurations encountered during 
the motion. 
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3. 1 . Qualitative physical picture 

Before proceeding, it is useful to present a simple qualitative picture which illustrates 
the idea of how entanglement can be established through classical motion even in 
noisy environments. We consider the following simple process: at time t = 0, two 
spins are far way from each other and the spin-spin interaction is very weak compared 
to the local magnetic fields. Thus the system Hamiltonian can be approximated as 
H M (0) ~ B(0)(aP + af ] ). Due to the presence of the thermal bath, the molecule is 
driven into its thermal equilibrium state with the corresponding ground state population 
denoted as p . If the local energy level splitting cj (0) = 25(0) is comparable or larger 
than the thermal energy scale, huo(0) > &bT, the ground state population p at this 
distant configuration will be relatively larger than the corresponding populations of 
the other energy levels. Nevertheless entanglement does not appear since the ground 
state is separable (|eo(0)) ~ |H))- Let us further assume that the spins come closer at 
such a speed that the coherent evolution due to the system Hamiltonian is adiabatic. 
While the spins approach each other, their interaction becomes stronger and the energy 
eigenstates adiabatically change into the ones of H M {t). In particular the ground state 
|eo(0)) — > |eo(t)) becomes more and more entangled as the coupling strength between 
spins increases, as indicated by the fact that its concurrence C(\eo(t))) = J/£ is an 
increasing function of the ratio J(t)/B(t). Following this reasoning, it is possible that 
entanglement appears as a transient phenomenon as long as the adiabatic change of Hm 
is still fast compared to the thermalization time scale 7 -1 with 7 = noj(2Nu,p + 1) for 
the chosen simple Ohmic spectral density. 



0.85 




50 100 150 200 250 300 50 100 150 200 250 300 
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Figure 4. Entanglement C (left) and the ground state population po (right) vs. t for 
a single conformational change with B{t) and J{t) as in Eq.|l8|l9]) and B(0) = 1.2, 
J(0) = 0.1, 5 = 0.01, to = 100, with different system-bath coupling strength (from top 
to bottom: n — 0, 0.002, 0.004, 0.006, 0.01). Here and in the following, the reservoir 
inverse temperature is j3 = 1. 

To make the argument more concrete, in Fig. (|4]) we illustrate the above qualitative 
picture through a simple process, in which 



B(t) 



B(0)-5-t 0<t<t 
B(t ) t < t 
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J<0) + *.| 0< ( < (0 

' J(t ) t„<t K ' 

Considering a reservoir with inverse temperature /3 — 1, the initial ground state 
population po is about 85%. The change of the system Hamiltonian is characterized by 
the small step value 5 = 0.01 which ensures the adiabatic condition for closed systems. 
It can be seen that the molecule can be kicked out of its separable thermal equilibrium 
state through classical motion and thus entanglement grows before being washed out 
by thermalization. The dissipative environment now plays its usual destructive role. 
Moreover, if the thermalization rate is fast compared to the adiabatic change of the 
system Hamiltonian, entanglement is greatly suppressed. 

3.2. Persistent recurrence of fresh entanglement 

To counter-act the detrimental effect of the environment and to sustain the generation 
of entanglement as a non-transient phenomenon, a persistent supply of free energy 
is needed. This is an organized form of energy which, as in our mechanism, can be 
introduced by conformational changes. To illustrate our idea, we first consider a two- 
spin molecule subjected to a deterministic and periodic oscillatory motion (the strict 



periodicity will be relaxed in section 3.4). The position of the two spins are given by 



2ixt 

xJt) = xJO) + (-l) Q a(cos 1) a = 1, 2 (20) 

T 

where x a (0) are the initial positions, a is the amplitude of oscillation, and r is the 



oscillation period. As a remark, we refer to the next subsection |3.4| for more general 
forms of classical motion involving stochasticity. The coupling between two spins is 
described by a dipole-dipole Ising interaction 

J(t) = J /d 3 (t) with d{t) = \xt(t) -x 2 (t)\. (21) 

while the local magnetic field profile, dependent on the spin position, is chosen as 

B a (t) =B - B ie - x «W a = B - 5 ie - d2 «/ 4CT , (22) 

The choice of x a (t), B a (t) and J(t) is quite arbitrary. As is explained in the following 
analysis, the essential condition for the proper function of the mechanism is that the 
dominant contribution to #m(£) is due to the background field when the spins are far 
away, and due to their interaction when they are closest. This is a natural assumption, 
since the interaction generally increases when the distance is reduced. 

We assume that the two spins start from the distant configuration and are initially in 
the corresponding thermal equilibrium state. Notice that this particular initial condition 
is irrelevant for the system asymptotic dynamics that we are interested in. At time 
t = the spins start to oscillate, i.e. they periodically come close to each other and 
separate again. The state of the molecule is driven out of thermal equilibrium, and 
entanglement is generated. When the spins move back to their distant configuration, 
the thermalization tends to increase the population of the ground state (due to the 
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increasing energy gap between the ground state and the first excited state), and thus 
effectively resets the system to a favorable condition. Therefore, although environmental 
noise generally destroys entanglement, it now represents a crucial factor to enable the 
persistent and cyclic generation of fresh entanglement. For simplicity, we call the cycle 
during the time interval [nr, (n + l)r] with n — > oo as the asymptotic cycle. It can be 
seen from Fig. ^ that entanglement indeed persistently recurs during the oscillations 
even if the thermal bath is so hot that no static entanglement can survive. 




20 25 30 35 40 20 25 30 35 40 



d d 



Figure 5. Left: Entanglement C vs. the molecular configuration d (the distance 
between the two spins) for the bosonic thermal bath with (3 = 1 and k = 0.01. The dot- 
dashed, blue curve indicates the first period; the solid, red curve marks the asymptotic 
cycle; while the green, dotted curve refers to the thermal equilibrium case. The arrows 
along the curves indicate how the entanglement evolves during the classical oscillations. 
The oscillation parameters are X\{Q) = —2:2(0) = —20, a = 5, r = 100, and Bo = 1.3, 
B 1 = 2.4, ct = 120, J = 10 4 . Right: Entanglement during the asymptotic cycle for 
different system-bath coupling strength (from top to bottom k = 0.005,0.01,0.02). 

In Fig. ^ we plot the population of the instantaneous ground state of the system 
Hamiltonian H.M{t) during the oscillation and observe how the periodic oscillations 
keep the molecule far away from its thermal equilibrium (as compared with the green 
curve). In particular, when the spins get closer and the entanglement of the ground 
state increases, the ground state occupation is noticeably higher than the thermal one, 
and this actually constitutes the key of the mechanism. Let us point out once more that 
thermalization tends to drive the molecule close to its instantaneous equilibrium state, 
but the continuous change of its shape prevents a full thermalization. The intrinsic reset 
effect of environmental noise is exploited to increase the ground state population when 
two spins reach the distant configuration. One can conclude that, even in presence of 
a very hot environment, there is still a chance to see entanglement in an oscillating 
molecule after an arbitrarily long time, see Fig. ([5]). 

3.3. Beat dephasing via thermalization 

The constructive role played by the environmental noise can be emphasized by including 
an additional source of decoherence to further destroy entanglement. Let us consider a 
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Figure 6. Population p of the instantaneous ground state of Hm(£) vs - the molecular 
configuration d. the oscillation parameters and the color code for the curves are the 
same as in Fig. ([5]). Here, the arrows along the curves indicate how the ground state 
population evolves during the classical oscillations. 



"phenomenological" pure dephasing process acting on the two spins via the term 

2 

V p p{t) = lp Y}^ P{t)af - p{t)l (23) 
i=i 

with 7 P denoting the dephasing rate. As it is well known, pure dephasing destroys all the 
coherence and always opposes the generation of entanglement. Even for the oscillating 
molecule, the effect of pure dephasing alone is to remove the off-diagonal density matrix 
elements of the two spin state, and thereby destroy entanglement in the long run. This 
is evident in Fig. ([7]), in which entanglement can be found only in the first oscillation 
period and certainly not in the asymptotic cycle (i.e. for t — >■ oo). What happens in the 
presence of a hot thermal bath? One would expect that entanglement were destroyed 
faster due to this extra decoherence source, and this is indeed true as long as only the 
first oscillation period is taken into account, see Fig. ([7]). However the constructive effect 
of thermalization becomes evident in the successive cycles in which it partly counteracts 
the pure dephasing through its built-in reset mechanism of re-pumping the ground state 
population (as explained in the above subsection). 

3.4- Stochastic classical motion 

So far we have considered a deterministic and purely periodic classical motion. In the 
context of biological systems, a more realistic approach naturally includes a certain 
degree of stochasticity. Different scenarios are conceivable depending on the biological 
situations that one intends to model. 

3.4-1- Fluctuating configurations. The allosteric process depicted in Fig. [T] is usually 
triggered by the interaction with a second chemical compound: when the molecule is 
in the open configuration, a chemical docks on the molecule and causes its allosteric 
contraction. The molecule reaches its closed configuration and stays there until the 
chemical is released. When the chemical undocks the molecule opens EZ]- 
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Figure 7. Entanglement C (left) and ground state population po (right) vs. the 
molecular configuration d. The dephasing rate is "f p = 0.02. The system-bath coupling 
strength is k = (purple curves) or k = 0.01 (red curves). The oscillation parameters 
are the same as Fig. ([5| left. 

The waiting times in the open and closed configuration can be thought as two 
stochastic variables which are, in principle, different. Supposing normal (Gaussian), 
but independent distributions for both processes, one obtains trajectories like the 
one depicted in Fig. (j8|, where the linear change in the spin distance during the 
conformational process is an arbitrary, but nonessential choice. 
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Figure 8. Example of a trajectory for waiting times distributed according to normal 
distributions (neglecting their tails for negative values). Mean d and standard deviation 
a in the open (closed) configuration are d = 14.5 , a = 6 (d = 7.5 , a = 5) respectively. 

Owing to our understanding obtained in the deterministic case, the intuition 
correctly predicts the limiting cases of long and short waiting times: when the chemical 
frequently docks or undocks from the molecule, the system is continuously kept out 
of equilibrium and its dynamics closely resembles the deterministic case. On the other 
hand, for long waiting times (i.e. long compared to the thermalization rate) equilibrium 
is reached in the extreme configurations and during the contraction of the molecule the 
single-shot behavior is reproduced. The creation of entanglement is here only a transient 
effect. 
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3.4-2. Potential landscape and thermal vibrations. A somewhat more interesting 
dynamics is obtained in the following scenario: the molecule has two configurations 
which minimize the potential energy associated to its classical shape. In our model 
it is immediate to identify the distance between the spins with an effective collective 
position coordinate for the macromolecule: varying such coordinate we move from the 
low-energy closed configuration to the low-energy open configuration via an energy 
landscape that contains barriers. For the sake of simplicity we model such energy 
surface as presenting only one central barrier and we add two steep increasing gradients 
to forbid configurations which are too extended or contracted, see Fig. (T9T) . 
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Figure 9. Energy landscape V(d) of molecular configurations with respect to the spin 
distance d. For consistency, we also set the temperature of this environment to (3 = 1 
and choose the barrier height accordingly. 



Once a configuration is chosen (as an arbitrary initial condition), the molecule starts 
exploring the other possible configurations due to the thermal energy supplied by the 
surrounding environment (notice that it is a totally different and unrelated environment 
with respect to the bosonic heat bath). It means that the molecule can randomly 
contract or stretch, as long as the energy of the new configuration is within thermal 
fluctuations. Using Metropolis algorithm and the potential landscape of Fig. ([9]), we 
obtain trajectories as the one shown in Fig. (10). 

In this situation, one cannot describe the molecular motion as composed by 
distinct opening and closing processes. In fact, the molecule stretches and contracts 
its configuration only partially and in a stochastic way. Nevertheless, the mechanisms 
presented in subsection are still at work in this scenario. Due to the increasing energy 
gap, the ground state population raises when the distance between the spins increases, 
with an "inertial" reaction time determined by the thermalization rate. Whenever the 
system undergoes a relatively fast contraction after having spent enough time in an 
extended configuration, entanglement is generated, as clearly appears in Fig. (11). 




Figure 10. Left: Example of a trajectory for the potential landscape depicted in 
Fig. ([9]) once thermal fluctuations are taken into account. Right: Magnification of a 
part of it. 
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Figure 11. Ground state (upper, blue) and concurrence (lower, red) in presence of 
the thermally driven fluctuations of the molecule's configuration. 



3.5. Non-Markovian case 

The microscopic derivation of the Lindblad type master equation previously introduced 
in this section relies on several approximations (for a more detailed derivation see [2Tj ) . 
These approximations refer to assumptions regarding the time scale of the thermal 
bath's dynamics and the strength of its interaction with the system. In particular, such 
interaction is supposed to be weak and the bath's relaxation time very short compared to 
the time scale of the system dynamics. For many biological processes, these conditions 
may not be satisfied. Nevertheless they can be considered as useful working hypothesis 
to obtain analytical expressions (suitable for numerical integration), and to illustrate 
our main idea. 

To demonstrate that such simplifying assumptions are indeed not essential for the 
validity of our arguments, we now extend our analysis to situations with a not-so- 
weak system-bath interaction and a non-Markovian reduced dynamics. We rephrase our 
problem within the path integral formulation of quantum mechanics: In this formulation, 
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the actual evolution of a quantum state is determined by the coherent sum of all the 
possible evolutions weighted via a phase proportional to the corresponding (classical) 
action. The presence of an environment in contact with the system leads to an additional 
weighting factor, the so-called influence functional [281 121] • The influence functional 
depends on the system's possible evolution and on the properties of the bath, and 
in general is non-local in time. This formalism is particularly convenient in case of 
environments composed of non-interacting harmonic oscillators, where one can obtain 
an analytical expression for the influence functional |29j . 

The part of the summation involving the system's degrees of freedom has to be 
performed numerically, and a particularly favourable prescription was introduced by 
Makri & Makarov under the acronym of QUAPI (QUasi Adiabatic propagator Path 
Integral) [3U1 EI]. Their method allows one to include a finite time memory via a 
truncated influence functional which takes into account only a few of the last propagation 
steps. The memory time is determined by how many time steps are included and by 
how long each of them lasts: T mem = Ak At. In principle, the longer T mem the further 
one can relax the Markovian assumption. In practice, however, the computational effort 
scales exponentially with Ak and the error in the single-step propagation quadratically 
in At. To overcome these limitations, each simulation has to be performed with 
different parameters and one has to perform two extrapolation procedures: the first 
keeps constant T mem while At — > to obtain precise estimate with finite memory, the 
second increases T mem — > oo to extrapolate the infinite memory time limit. 

These technical procedures are extensively presented in [32] and are beyond 
the aim of this work: Here we do not want to precisely quantify the amount of 
entanglement generated under a specific system-environment interaction, but our goal is 
to demonstrate that the mechanism we introduced is indeed robust and is not confined 
to the Born-Markov regime. However, we have selectively performed the extrapolation 
procedures for several points along the curves: The agreement is always good (on the % 
level) and consistent with the values of entanglement found before the extrapolations. 
The extent of the memory effects is determined by the so called bath response function 
BRF: To properly take into account the non-Markovianity of the environment, the 
memory time T mem has to be larger than the bath correlation time given by the width 
of the BRF, see Fig. (12). Fig. (13) shows how the generation of entanglement persists 
even for baths with finite-time correlations (determined by a small value for the cut- 
off frequency in the spectral density J(u), giving rise to finite memory effects) and 
moderately strong system-environment couplings. 



4. Spin gas model 

The actual environment for biological systems is expected to be very specific and 
complex, and has not yet been fully characterized. In the above section we investigated 
the case of thermal baths composed of a large number of harmonic oscillators, but this 
is most probably still a simplistic (even if reasonable and common) choice. In this 
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Figure 12. Bath response function for Markovian (us c — 40, top) and non-Markovian 
(lj c = 4, bottom) regimes. The solid, red curve is the real part of the BRF, the dashed, 
green curve is its imaginary part, while the vertical, black lines mark the memory time 



1.5 considered in Fig. (13) 
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Figure 13. Entanglement C (left) and ground state population po (right) vs. 
molecular configuration d for a non-Markovian environment. The dot-dashed, blue 
curve indicates the first period; the solid, red curve marks the asymptotic cycle; 
while the green, dotted curve refers to the QUAPI equilibrium state for each specific 
configuration. Top: uj c = 4, n = 0.01, T mem = 1.5, r = 100, At = 5; Bottom: 
lo c = 4, k = 0.05, T mem — 1.5, t = 75, At = 5; while the other parameters are as 
in Fig. ([5} left. 



perspective, it is worth to see whether the main features that we observe for bosonic 
heat baths are preserved also for a different kind of environment, namely the spin gas 
model [321 [Ml ES], which is particularly suggestive in a biological context. Here the 
environment is composed by particles which follow random classical trajectories and 
collide with our two spins. During the collisions they stochastically interact with the 
spins giving rise to an effective decoherence process. Under certain assumptions [35J, 
at a coarse time level that includes both the coherent Hamiltonian dynamics and the 
random collisions, it is possible to derive a phenomenological description in terms of a 
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Lindblad-type quantum master equation characterized by the following dissipator: 

2 

Vp{t)= 1 Y}^ ) P{t) + C^ ) p{t)l (24) 

a=l 

with 

K g'p = s[2a y + 'pa_ ' — pa_ cr\ — a_ 'a\ ' p\ (25) 

r( a ) /i \ To (<*) (ot) M M (a) (a) i /nr\ 

£ y d ' p = (1 — s)[2a_ pa + — pa + <r_ ' — a y + 'a_ 'p\, (26) 

where a± = (ai^ ± io r y a ' ) )/2 are the ladder operators for a two-level system. The de- 
coherent channels Cg (p), £-d\p) represent the energy gain and loss during the process: 
7 is the rate with which ((1 + cr 2 )/2) decays, and s is its equilibrium value Without 
loss of generality we assume that 7 > and < s < 1/2. We also note that correlations 
in the spin gas can be included in the analysis [37] . 



4-1- Static entanglement in the steady state 



Similar to the case of the bosonic heat baths, we first consider the steady state 
q = p{t — > 00) for each molecular configuration, i.e. at fixed spin-spin interaction 
strength J and background field B, to verify that it is a separable state. To this end, 
we impose 

-i[H M (t),p]+T>p = 0. (27) 

It is straightforward to check that the steady state for each configuration d is unique 
and that its non zero density matrix elements are given by 

6m = [J 2 + 46 2 s 2 ]/4(J 2 + 6 2 ) 
gn ^[J 2 + Ab 2 s(l-s)}/A(J 2 + b 2 ) 
te = [J 2 + 46Ml- S )]/4(J 2 + 6 2 ) 
£33 = [</ 2 + 4& 2 (l- S ) 2 ]/4(J 2 + 6 2 ) ' 
g 03 = -J(2B - i 7 )(l - 2s)/2(J 2 + b 2 ) 
£30 = -J(2B + i 7 )(l - 2s)/2(J 2 + b 2 ) 

[AB 2 + 7 2 ) 1 / 2 . The static entanglement of the steady state quantified by 



where b 

concurrence is calculated as C(g) 



2max{0, |^ 3| - (^iiP22) 1/2 }, resulting in 



C{g) = 2max{0. 



2J6(l-2s)-[J 2 + 46 2 (l-s)s] 



}• 



(29) 



4(J 2 + 6 2 ) 

Thus, no static entanglement is possible if the parameter s exceeds the critical value of 



max{0, ^[1 + J/6 - (1 + 2J 2 /6 2 ) 1 / 2 ]} 



(30) 



plotted in Fig. (14). 
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Figure 14. The critical value of s (above which no static entanglement exists in the 
steady state) for different values of B and J with 7 = 0.025. 



4-2. Persistent dynamic entanglement 

The physical picture underlying the spin gas model is very different from the collection of 
harmonic oscillators constituting the bosonic thermal bath and these two environments 
destroy entanglement according to different mechanisms. Nevertheless, they share 
the essential features which allow the cyclic generation of fresh entanglement via the 
constructive role played by the environment noise itself. Even if in any possible static 
configuration the steady states for both models and high enough temperature (bosonic 
heat baths) or level of noise (the spin gas environment with a large value of s) does not 
exhibit entanglement, it persistently appears when classical oscillations are taken into 
account. 

This can be explained since the spin gas environment also contains a built-in 
reset mechanism via spin flip: For the molecular configuration in which the two spins 
are spatially separated, the environment effectively lowers the system's entropy and 
increases the ground state population. Thus, under these conditions and in accordance 
with what has been observed for thermal baths, dynamic entanglement is generated in 



the long time limit, Fig. (15) 



As in the case of the bosonic heat baths, here the constructive role of the 
environment noise can also be further illustrated by adding pure dephasing. One can 



see from Fig. (16) that the noise in the spin gas model can counteract the effect of 



pure dephasing and sustain the persistent dynamic entanglement. The green curve is 



here obtained from an expression analogous to Eq.(28), and generalized to include the 
additional dephasing. 
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Figure 15. Entanglement C (left) and ground state population p (right) vs. the 
molecular configuration d for the spin gas environment. The dot-dashed, blue curve 
indicates the first period; the solid, red curve marks the asymptotic cycle; while 
the green, dotted curve refers to the steady state for a specific configuration. The 
decoherence rates are s = 0.16 and 7 = 0.025; the oscillation parameters are the same 
as in Fig. ^ left. 
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Figure 16. Entanglement C (left) and ground state population p (right) vs. the 
molecular configuration d for the spin gas environment with additional pure dephasing. 
The dot-dashed, blue curve indicates the first period; the solid, red curve marks the 
asymptotic cycle; while the green, dotted curve refers to the steady state for a specific 
configuration. The decoherence rates are s — 0.16, 7 = 0.025, and j p — 0.01; the 
oscillation parameters are the same as in Fig. jjj) left. 

5. Concluding remarks 

In this work, we have demonstrated how entanglement can persist as a recurrent property 
even in systems exposed to a hot and noisy environment. To achieve this task we have 
exploited the presence of an internal driving force which keeps the system out of the 
corresponding equilibrium state. In biologically inspired context, such driving force can 
be identified with conformational changes, including the stretching of macromolecules 
or their folding/unfolding processes. Here, we describe these conformational changes as 
purely classical, but capable of influencing the dynamics of localized quantum degrees of 
freedom and to prevent complete thermalization. The non-equilibrium interplay between 
classical motion and open quantum system dynamics gives rise to an effective reset 
mechanism which promotes the generation of fresh entanglement even for arbitrary long 
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times. The present model serves two purposes: On one hand, it provides a concrete 
route where to search for signatures of entanglement in biomolecular systems operating 
at room temperature. On the other hand, it indicates new directions to design molecular 
scale machines able to generate non-trivial quantum states even in noisy environments. 
Finally, let us remark that an experimental simulation of our model is achievable with 
state-of-the-art technology in setups like ion traps or nano-mechanical systems. 
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